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Abstract 

For the first time, we study the evolution of the stellar mass-size relation for star- forming galaxies 
from z A io z 7 from Hubble- WFC3/IR camera observations of the HUDF and Early Release 
Science (ERS) field. The sizes are measured by determining the best fit model to galaxy images 
in the rest-frame 2100 A with the stellar masses estimated from SED fitting to rest-frame optical 
(from Spitzer/IRAC) and UV fiuxes. We show that the stellar mass-size relation of Lyman-break 
galaxies (LBGs) persists, at least to z ~ 5, and the median size of LBGs at a given stellar mass 
increases towards lower redshifts. For galaxies with stellar masses of 9.5 < log(M*/M0) < 10.4 sizes 
evolve as (1 + ^^-i.2Q±q.ii ^ This evolution is very similar for galaxies with lower stellar masses of 
8.6 < log(M^/M0) < 9.5 which is rg cx (1 + ^)-i i8±o.io^ -^^ agreement with simple theoretical galaxy 
formation models at high z. Our results are consistent with previous measurements of the LBGs 
mass-size relation at lower redshifts ~ 1 — 3). 

Subject headings: galaxies: evolution - galaxies: high redshift - galaxies: structure 



1. INTRODUCTION 

The size of a galaxy is a fundamental and important 
parameter to measure. Over the past decade, observa- 
tions have revealed that sizes of galaxies at a given stellar 
mass were smaller at higher redshifts and change signifi- 
cantly with redshift. It has been shown that the sizes of 
galaxies correlate with their stellar masses and that this 
correlation exists at least up to 2; ^ 3 (g-R-, iFranx et al 



20081: IWilliams et al.ll2Mol : lMosleh et al.ll2011l : [Law et al 



20121 ) 



There are many proposed scenarios to explain the 
physical processes of galaxy assembly that plausibly re- 
produce the observable stellar mass and size of galax- 
ies at di fferent redshifts ( e-g-, galaxy minor or major 
mergers fKhoch far fc^iffl [20061 [20091 : iBell et al.l [200I 
iNaab et"a l. 2009) or g:as accretion in outer regions an d 
star formation (Dekel et al.l l2009l :lElmegree n et al.l l2008l. 
Accurate measurements of both stellar masses and sizes 
of galaxies over a wide redshift range are fundamentally 
important to constrain these galaxy formation models. 

To extend the mass-size relation of galaxies to the 
highest redshifts, we exploit the Lyman-break galax- 
ies (LBGs) which are star forming galaxies with strong 
rest-frame UV emission and could be selected by pho- 
tometric dropout tech niques (e.g. Steidel et al. 20031 ; 
lAdelberger et al.l l2004f ) . These galaxies c an be identi- 
fied out to very high redshifts (e.g., z ~ 8)C0 esch et al.l 
I2OI2I : lYan et aU |2011[ ) and thus provide insight into the 
early evolution of the mass-size relation. 

Morphological studies of LBGs (2: ~ 2 — 6) in rest- 
frame UV have shown that these galaxies are mostly 
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compact sources however, multiple core LBGs have also 
been found (e.g. [R avindranat h et al.l [2006; L aw et all 
120071 : IConselice fc Arnold 2009). Analysing their size 
and structure could help to interpret the dominant mech- 
anism for galaxy growth. The new Wide Field Camera 
3 (WFC3) on board the Hubble Space Telescope (HST) 
can provide sizes of high-z galaxies in longer rest-frame 
wavelengths than Advanced Camera for Survey (ACS). 

Size studies of galaxies at redshifts z > 4 usin g 
profile fitting techniques are rare (Oesch etHI 120103) . 
Here for the first time, we investigate the mass-size 
relation of dropout galaxies up to the very early stages 
of galaxy formation, using the advantages of wide-area 
HST surveys and high spatial resolution of WFC3. 
We measure the sizes of LBGs at approximately the 
same rest-frame wavelength at different redshifts, 
minimizing the effects of morphological K-correction. 
We utilize observations of both Hubble Ultra Deep 
Field (HUDF) and the deep wide-area Early Release 
Science (ERS) field to study the mass-size relation of 
the largest sample of LBGs at 2: ~ 4 — 7 so far. The 
cosmological parameters adopted throughout this pa- 
per are = 0.3, = 0.7 and Hq = 70 km Mpc~^. 

2. DESCRIPTION OF DATA 

The sample of dropout sources at z ~ 4 — 7 used 
here is based on recent HST-WFC3/IR observations 
over the HUDF and ERS field ( see Bouw ens et al.ll2011l : 
lOesch et al.ll2010bl : [Gonzalez et al. 2010) . Ou r samp le of 
dropout galaxies is tak en from [Bouwens et alj (|2011 ) and 
[Gonzalez et al.l (|2011l ). Their candidates are selected by 
utilizing color criteria (Lyman-Break technique) for se- 
lecting B, V, i and z dropouts simi l ar to those used in 
lOesch et al.[ (|2010br ): [Bouwens et all (j2 0071 [2011 1 ) . Their 
sample consists of 679 objects at z ^ 4 — 6 in ERS field 
and 345 dropouts at z 4 - 7 in HUDF. 

For all candidate galaxies in this paper we us e 
the stellar masses measured by .Gonzalez et al.[ (|2011[ ) . 
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Fig. 1. — From left to right, H-band postage stamps (6^^ x 6^^), 
best-fit models from GALFIT, residual images and mask maps are 
shown for a z-dropout candidate (top row) and an i-dropout can- 
didate (bottom row). The quality of the fits can be seen from the 
residual images. 

These authors used FAST (jKriek et alJ I2QQ9D to 
fit template galaxy SEDs to photometry in ACS, 
WFC3/IR and IRAC [3.6 ] and [4.5] filters. They used 
iBruzual fc Charlotl (|2QQ3f ) stellar population evolution 
models w ith const a nt sta r formation histories and as- 
sumed a iSalpeterl ()1955f ) initial mass function (IMF; 
0.1 - IOOMq) and 0.2^0 solar metallicity. Reliable de- 
blended IRAC fluxes and photometry in different bands 
(with different PSFs) were obtained by using a source- 
fitting algorithm d escrib ed in iLabbe et al.l (j2006|) (see 
also Gonzalez et an[2QlQf ). 

The photometric deblending for the IRAC imaging 
does not properly work for objects in crowded regions due 
to the large PSF of IRAC bands; due to this, the number 
of sources with re liable m ass estimates is reduced to 60% 
(see also Gonzalez et al.l '2011). The stellar masses are 
corrected to a KrouDa (2001) IMF by a reduction of 0.2 
dex (see 'Marchesini et al..>2009.) . These are consistent 
with the SDSS masses. 

The near-IR images used in this study are taken 
from the full two-year WFC3 /IR H UDF (Yio5,Ji25, 
ij^i^n: (|Bouwens et al.l 120111: lOesch et al.l l2010bD ) and 
ERS (Yogs, ^125, i^ieo; (iB ouwens et~alll2011f )) data. The 
data set covers ^ 4.7 arcmin^ in the HUDF and ~ 40 
arcmin^ in the ERS (in GOODS-South) field with a 
pixel scale of 0.06'^ and a point spread function (PSF) 
FWHM of - 0.17'' in the H-band. We also use very 
deep images of t he HUDF field obtained with the ACS 
(jBeckwith et al.l [2006) and deep GOODS ACS/ WFC 
data over the GOODS field C Giavalisco et"lal[200l with 
a pixel scale of 0.03'' and a PSF FWHM of ^0.10". The 
reduction of these data is described in iBouwens et al.l 
(j2011 ) and references there in. 

As explained in the next section, we use sizes of 
objects in the band closest to the rest-frame 2 100 A. 
Therefore, in addition to reliable stellar mass cut, we 
restrict our size measurements to objects with sufficient 
S/N in their size measurement band. Hence, we use 
HiGo < 28.5 in HUDF and Hiqq < 26.2 AB mag in 
ERS field (corresponding to lOa in 0.5" apertures) 
for i-dropouts and z-dropouts, J125 < 28.5 in HUDF 
and J125 < 26.5 AB mag in ERS for V-dropouts and 
Yios < 28.3 in HUDF and lo98 < 26.4 AB mag in ERS 



for B-dropouts. Our simulations (see Section 3) show 
that systematic uncertainties on size estimates are small 
up to these magnitudes. Moreover, objects with poor 
S/N (< 5cr) are not included in our analysis. Thus, 
the final sample consists of 156 B-dropouts {z ^ 4), 
45 V-dropouts {z ~ 5), 13 i-dropouts {z ^ 6) and 4 
z-dropouts {z ^ 7). 

3. SIZES 

Ou r size measurements are performed in a similar way 
as in "Mosleh et al.l ()201 1"). PSF-convolved Sersic models 
are fit to all g alax ies using the GALFIT version 3 package 
(jPeng et al.ll20Toh . We measure sizes in filters described 
in the previous section. 

GALFIT requires an accurate point spread function 
(PSF) of the image to convolve the Sersic profile mod- 
els in order to find the best (x^ minimized) fit. For 
the HUDF W FC3/IR ima ges we use a PSF created us- 
ing TinyTim (jKristI 1 19951). Although, we prefer using 
empirically-determined PSFs (from stars), in the HUDF 
there are not enough stars to create a sufficiently low- 
noise PSF for deconvolution. We note that our test 
shows that sizes measured using non-saturated star as 
a PSF give results which are consistent with those based 
on TinyTim created PSFs. For the ACS images in the 
HUDF we use non-saturated stars in the field to make a 
median-stacked PSF. In the ERS field, we use a median 
stacked-PSF using non-saturated stars in the field for all 
optical and near-IR bands (ACS and WFC3/IR). 

GALFIT measures the half-light radius along the semi- 
major axis, a, and the axis ratio, 6/a, of each galaxy. We 
determine the circularized effective radius of the galax- 
ies (re = a^hjd) to remove the effects of ellipticity. We 
fix the Sersic index to 2.5 where GALFIT fails to return 
reliable measurements; i.e., for objects with large uncer- 
tainties in their output parameters. Size determination 
for all objects while fixing the Sersic index to 2.5 re- 
veals that this will not introduce systematics in our size 
study. Fixing Sersic index to 1.5 or 3 gives similar size 
estimates as for n = 2.5. In Figure [T] we show from left 
to right, postage stamps in the H-band, best fit models 
from GALFIT, residual images and mask maps for a z- 
dropout (top row) and an i-dropout (bottom row) in the 
HUDF. The low level residuals shows that our best fit 
models closely match the observed galaxy images. 

We first test the accuracy of our measurements by com- 
paring the effective radii, rg, measured at different wave- 
lengths. A comparison of the sizes of our dropout can- 
didates is shown in Figure [2l Each dropout sample is 
represented by a different color: B-dropouts are black, 
V-dropouts are green and i-dropouts are blue. We note 
that ACS has a much smaller PSF and finer pixel scale, 
and it is known to produce highly reliable size measure- 
ments from previous studies. The correlation between 
estimated sizes in different bands, is an indicator of the 
size measurements reliability. Since color gradients and 
clumpy star formation may introduce additional scat- 
ter between bands, these comparisons essentially show 
an "upper limit" to the intrinsic uncertainties. 

In addition, we have performed realistic simula- 
tions on both the H UDF and ERS fields analogous to 
'Mosleh et aJj (|2011 ^. by adding model galaxies to the 
images and measuring their parameters in the same way 



Evolution of the mass-size relation From z = 1 — 7 
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Fig. 2. — Comparison of sizes of dropout candidates measured in different filters with SjN > 10 in both the HUDF and ERS. Each 
color represents a different sample of dropout candidates. The correlation between estimated sizes in different bands is an indication of the 
reliability of our size measurements. The small offsets between sizes in different bands maybe due to color gradients in the galaxies. The 
objects with large offsets from the one-to-one relation have generally a clumpy structure or multiple cores. 



as for real objects. From the simulations, the relative 
difference between the recovered and input sizes versus 
magnitude shows that the systematic uncertainties on 
the recovered sizes are very small (< 10%) down to S/N 
ratio of 10 (corresponding to i^ieo < 28.5 in the HUDF 
and iJieo < 26.2 in the ERS). 

Finally, in order to minimize the effects of morpho- 
logical i(^-correction, we use sizes of each object in the 
WFC3/IR or ACS band that is closest to the rest-frame 
2 100 A. Therefore in this study the Yios band images 
are used for B dropouts, J125 images are used for V 
dropouts and i^ieo band images for i and z dropouts. 

4. RESULTS 

4.1. Size-mass Evolution 

The size and stellar mass estimates of our dropout can- 
didates allow us to investigate the mass-size relation and 
its evolution. Figure [3] shows the mass-size relation for 
four redshift bin from z = 4 to z = 7. There is a clear 
relation out to z ^ 5 and a hint at z ~ 6 for galaxies with 
log (M^/M©) > 8.6. In each panel, color points represent 
the distributions of candidates in the mass-size plane. 
Objects with S/N < 2(j in the IRAC [3.6] channel, are in- 
dicated by open symbols, and those with higher S/N (i.e. 
more reliable stellar mass) are shown as filled circles. The 
pink stars indicate the median sizes of galaxies in differ- 
ent stellar mass bins down to our estimated stellar mass 
limit. The dot-dashed line in each panel represents the 
mass-size relation for late- type galaxies (i.e., n < 2.5) in 
the local universe from Shen et al. (2003) and the dashed 
line represents the best-fit stellar mass-size rela t ion fo r 
UV-bright galaxies at z - 2 from iMosleh et all i\201l\ ). 
According to our simulations, sizes below the dotted lines 
(corresponding to apparent radii of 0.05'^) have large un- 
certainties and can be conservatively thought of as upper 
limits. 

The stellar mass-size relation for z ~ 4 and z ~ 5 is 
characterized by rg oc with a = 0.14 ± 0.06 and 
0.17 ± 0.07, respectively. These are significantly differ- 
ent from zero. We note that the correlation coefficients 
are 0.22 at z ~ 4 and 0.37 at z ~ 5. At z ^ 6 there 



might be a hint of a mass-size relation, but there are too 
few galaxies in our sample at this redshift to either con- 
firm or rule out a relation. We note that the size and 
stellar mass of these galaxies need not to be correlated 
as the sizes are measured in the rest-frame UV and the 
stellar masses in optical rest-frame. Nonetheless, the re- 
sults show the persistence of the mass-size relation for 
star- forming galaxies up to very high redshifts. We also 
note that the fraction of contaminants in our dropout 
samples is r elatively small (< 11%) (see iBouwens et al.l 
([2OIII . [2OO7I ) for more details) and hence unlikely biased 
the relations. 

The black diamond in the lower left panel of Figure 
[3] shows the position in the mass-size relati on of an i- 
dropo ut galaxy {SBM03#1 in HUDF) from lEvles et all 
(|2005f ). This object has a robust spectroscopic redshift 
{z = 5.83) and a stellar mass of a few times lO^^M©. 
The size estimate for this object is consistent with the 
median size estimate of i-dropouts in our higher mass 
bin. 

4.2. Size Evolution at fixed mass 

Comparing the median sizes of LBGs for a given stel- 
lar mass at different epochs, illustrates the evolution of 
galaxies effective radii at fixed stellar masses with red- 
shift. This size evolution is shown in Figure IH Galaxies 
are split into two stellar mass bins: 9.5 < log (M^/M©) < 
10.4 (shown in the left panels) and 8.6 < log(M*/M0) < 
9.5 (shown in the right panels). For the objects in the 
lower redshift range (i.e. z ~ 1 — 3), we used sizes mea- 
sured for UV-bright galaxies in the same stella r mas s 
range from our previous studies by Mosleh et a Lj (|2011l ). 
In order to measure the size evo lution consistently we 
use similar technique described in iNewman eFall ([2QT2h 
to normalize sizes of galaxies with re cx: M^-^^ to stellar 
mass of 10^- (left panels) and lO^M© (right panels). 
In the bottom panels, blue and green points represent 
normalized half-light radii of dropout objects as a func- 
tion of redshift. 

We fit a simple power law of the form (1 + z)~^ to the 
observed points (bottom panels of Figured]); the best 
fits are shown as solid lines. In the high stellar mass bin. 
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Fig. 3. — The mass-size relation for dropout galaxies in different redshift bins. The filled circles indicate sizes for galaxies with reliable 
stellar masses (> 2cr detection in IRAC[3.6]). The pink large stars are the median of sizes from GALFIT in different stellar mass bins with 
the errors calculated by bootstraping.The black diamond in the lower left panel is an i-dropout galaxy at z ~ 6 with a stellar mass of a 
few times VS^^M(^ from Evles et al. (2005). The dashed and dot-dashed lines show the mass-size relation for UV-bright galaxies at z ~ 2 
from Mosleh et al. ( 2011) and for late- type galaxies at z ~ from Shen et al. ( 2003), respectively. The dotted lines correspond to our size 
measurement limits. Typical errors of the points are shown in the lower right of each panel. This plot shows that mass-size relation holds 
out to 2; ~ 5, and possibly beyond for galaxies with log(M*/M0) > 8.6. Sizes decrease toward higher redshift. 



galaxy sizes are found to evolve as (1 + ^) i-20±0-ii^ and 
for galaxies with lower stellar masses the size evolves as 

(1^^)-1.18±0.10^ 

In the top panels of Figure [H blue and green soHd 
circles represent the median effective radii of dropouts 
from this study, and the solid squares are median sizes 
in lower redshift bins. The r ed poi nts are the mean 
galaxy sizes from lOesch et alJ (|2QlQa[ ) for two different 
luminosity ranges: (0.3 — V)L^ ^=3 in the left panel and 
(0.12 — 0.3)L* in the right panel. 

The filled yellow and open black stars in the top pan- 
els are sizes based on stacking galaxies from fixed stellar 
mass ranges for our two fields. For stacking, we used the 



central positions of objects determined by GALFIT and 
replaced contaminated pixels from neighbouring sources 
with sky background values. We then used GALFIT to 
measure the half-l i^ht radii of our fina l sta cked images. 
As shown by e.g., iHathi etldl (|2008a[ ) and lOesch et al.l 
(2010a), stacking can reproduce reliable average surface 
brightness profiles. Thus, the agreement between the me- 
dian points and results of stacking suggests that we are 
not systematically missing light in the extended wings 
of galaxies and that on average our size estimates are 
robust. 

5. SUMMARY &; DISCUSSION 
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Fig. 4. — The redshift evolution of the sizes of dropout galaxies in two stellar mass ranges: log(M*/M0) ~ 9.5 — 10.4 (left panels) and 
log(M*/M0) ~ 8.6 — 9.5 (right panels). The lower panels show the distribution of sizes for all sources in our analysis and the sizes of 
UV-bright galaxies at z ~ 1 — 3 within the same mass range from Mosleh et al. (2011). The best power law fits to these data points are 
shown as solid lines. The red points are from Oesch et al.. ( 2010a ) for a luminosity-selected sample. The black and yellow stars are the sizes 
derived by stacking galaxies in the same stellar mass range. The black diamonds in the top panels indicate median size of local late-type 
galaxies from Shen et al. (2003). 



For the first time, we have studied the stellar mass- 
size relation of LBGs out to 2; ~ 7 using ultra-deep 
WFC3/IR observations in the HUDF and ERS fields. 
We have shown that the mass-size relation of star- 
forming galaxies persists to very high redshift s, and 
that at fixed stellar mass, the sizes of galaxies increase 
significantly towards later cosmic time. The observed 
size growth of LBGs studied here - rg cx (1 + ^)-i-20±o.ii 
for galaxies with stellar mass of lO^'^-lO^^'^M© - is 
in agreement with previous stellar mass-size studies 
at z <. 3 (e.g;., iDahlen et alJl2QQ7l : iMosleh et ani2Qlll : 
iNagv et all 1201 ll : iLaw et al.ll2Q12l ). It is also consistent 
with the size evolution estimated by oth er studies based 



on l u minosity-selected samples f e.g.. iFerguson et al.l 
2QQ1: iBouwens et alJ [2QQ1 . [2QQ1 iHathi et all l2QQ8bl : 



Oesch et al.ll2QlQa[ ). 



The redshift dependence of the size evolutions is very 
similar for both high and low stellar masses. Therefore, 
the galaxy size evolution might be written as a separa- 
ble function of mass and redshift and this would be in 
agreement with simple models of galax y formation de- 
veloped for high redshift systems (see IWvithe fc Loebl 
[2011). However, our sample is not stellar mass com- 
plete; hence using deeper samples in future is needed for 
further investigations. 



In Figure [5] we compare our e stimated size-reds hift re- 
lation for LB Gs with those in iLaw et al.l (|2012f ) (blue 
triangles) and iDahlen et al.l (|2QQ7l ) (open squares). The 
sizes are normalized to a stellar mass of lO^^M©. The 
size estimates from different studies are consistent with 
the best fit found in this study (m = 1.20 ± 0.11, solid 
line). This suggests that LBGs may evolve into UV- 
bright galaxies at 2: ~ 1. At 2: ~ 6 these galaxies are 
extremely compact: rg ~ 0.8 kpc, at stellar mass of 
10^^ Mq. However, they grow by a factor of about 6 
to z ~ 0.85. The z ^ 6 galaxies have a mass size rela- 
tion close to those of the compact quiescent galaxies at 
z ~ 2; For example the normalized size (re/(M/10^^)^'^ 
) of a sample of quie scent galaxies at 2: ~ 2 studied by 
iSzomoru etHI (|2011f ) is - 0.6 kpc. 

The scarcity of observed LBGs between z ^ and 
2; < 1, complicates the interpretation of the evolution 
of their stellar mass-size relation to the present time. 
The median size measured for local late- type (i.e., 
n < 2.5) SDSS galaxies by iShen et al.l ^003) at the 
same stellar mass (black soHd diamond in Figures |4] and 
[5l corrected to t he rest- fram e UV using the analys is by 
iGil de Paz et al.l (j2007[ ) and lAzzollini etld] (j2009f ) and 
D. Szomoru (private communication)) is smaller than 
the size of z ^ 1 UV-bright galaxies. We note that 
the SDSS late-type sample is most likely a different 



6 



Mosleh et al. 



Redshift 

0.5 1 1.5 2 3 4 5 6 7 




0.0 0.2 0.4 0.6 0.8 1.0 

log(1 +z) 



Fig. 5. — Comparison of sizes of UV-bright galaxies measured by different authors. In order to have a consistent comparison, sizes are 
normahzed using Ve oc M^-^^ for a stellar mass of IO^^Mq. The solid line is the best fit from our analysis (re oc (1 + 2;)~^-20±o.ii^ r^j^^ 
dotted lines represent the original (not normalized) size measurement. The black solid diamond is the median size of late- type galaxies at 
z ~ for stellar masses around Mq from Shen et al] (|2003l ). with the measured dispersion shown as error bar. The sizes of UV-bright 
galaxies at fixed stellar mass increases rapidly towards later cosmic time. However, the sizes measured for local late-type galaxies are 
smaller than the sizes of UV-bright galaxies at 2; ~ 1. These samples most likely comprise intrinsically different galaxies. 



galaxy population than the UV-bright sources we 
study at z > 1, and is therefore not relevant for direct 
comparison. lOverzier et al.l (j2QlQf ) studied 30 local 
Lyman Break Analogs {z < 0.3), however this sample 
was selected to have a surface brightness limit and is 
therefore not characteristic of star-forming galaxies in 
the nearby universe. In addition, there is some evidence 
that the evolution of the stellar mass-size relation for 



star- forming galax ies is slower between z = 1 and 
z = (e.g.. Bar d en et al.l [2QQ5I ). This needs further 
investigation using homogeneously-selected sample in 
future. 
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